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We use a quasi Hamiltonian formalism to describe the dissipative dynamics of a circuit QED qubit 
that is affected by several fluctuating two level systems with a 1/ f noise power spectrum. The qubit- 
resonator interactions are described by the Jaynes Cummings model. We argue that the presence 
of pure dephasing noise in such a qubit-resonator system will also induce an energy relaxation 
mechanism via a fluctuating dipole coupling term. This random modulation of the coupling is 
seen to lead to rich physical behavior. For non-Markovian noise, the coupling can either worsen or 
alleviate decoherence depending on the initial conditions. The magnetization noise leads to behavior 
resembling the collapse and revival of Rabi oscillations. For a broad distribution of noise couplings, 
the frequency of these oscillations depends on the mean noise strength. We describe this behavior 
semi-analytically and find it to be independent of the number of fluctuators. This phenomenon 
could be used as an in situ probe of the noise characteristics. 



I. INTRODUCTION 

In the area of quantum information processing(QIP), 
cavity quantum electrodynamics(QED) has received 
much attention as it allows the inter conversion between 
matter based qubits and photonic qubits. There are sev- 
eral promising cavity QED schemes - ion traps 1 , neu- 
tral atoms 2 , Rydberg atoms^, quantum dots in photonic 
crystals 4 . However superconducting qubits coupled to 
transmission line resonators (or LC-circuit oscillators)P^l 
appear to be extremely promising as they offers two 
substantial advantages. Firstly, existing fabrication and 
lithographic techniques can be used to make a large 
number of scalable superconducting qubits, which act 
as macroscopic quantum objects. Secondly, the smaller 
mode volume of the resonators and the larger physical 
size of the qubits, allows strong coupling between the 
cavity mode and the qubit, enabling efficient manipula- 
tion of interactions at the single photon level. 

In cavity QED systems and particularly in the case of 
circuit QED schemes, there are three main mechanisms 
that contribute to the qubit 's decoherence. They are: 
decay into photon modes other than the cavity mode, re- 
laxation due to vacuum fluctuations and dissipation due 
to electrical noise. The first one is negligible for trans- 
mission line resonators largely due to their 1-D nature. 
The relaxation due to vacuum fluctuations can be re- 
duced by increasing the level detuning 9 . The third cat- 
egory can be further divided into extrinsic and intrinsic 
sources of noise. In case of SC flux/phase qubits, extrin- 
sic noise sources, such as ones due to the fluctuations in 
the external circuitry can, for example, be fixed by in- 
creasing the impedance of the current source providing 
the flux bias^ . On the other hand intrinsic noise of the 
1// typ d 11 * 12 1 , is often extremely difficult to deal with. 
However, it has been shown that two level fluctuators 
contributing to 1/ f noise can be minimize d by suitably 
engineering the Joshephson tunnel junctiori 13 * i 

While 1/f type magnetic f lux n oise was observed in 
SQUIDs quite some time ag d 15 * 16 * its origins remained 



somewhat unexplained. Recent interest in quant um co m- 
puting has however renewed interest in this fiel d 11 * 12 *. A 
number of models have been suggested to explain this 
phenomenon. Examples include a model proposed by 
Koch et al 17 where unpaired and non-interacting elec- 
trons randomly hop between traps with fixed but random 
spin orientations. The trap energies in their model have 
1/f distribution. A dangling bond modeP^ for 1/f noise 
has been proposed where electrons flip their spins due to 
their interaction with tunneling two-level-systems(TLS). 

Relatively recent experiments^ suggest that flux noise 
arises from unpaired surface spins which reside at the 
superconductor-insulator interface in thin- film SQUIDs. 
Choi et aP^have explained this experimentally estimated 
areal spin density in terms of metal induced gap states 
that arise due to the potential disorder at the metal- 
insulator interface. These same set of experiments have 
observed strong correlation between the inductance noise 
and the flux noise, suggesting that the surface spins 
are interacting 19 . This appearance of a type of a long 
range magnetic order is somewhat consistant with a 
model suggested by Faoro and Ioffe 2 ^ where the spins 




FIG. 1. Schematic showing a superconducting(SC) qubit, 
coupled to a resonator, under the influence of several fluctuat- 
ing two level systems present in the SQUID 's metal-insulator 
interface. 
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at the superconductor-insulator interface interact via the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism. 

Regardless of the exact microscopic details of the noise 
mechanism, there appears to be a consensus that the 
1/f type magnetization noise arises from the presence 
of TLSs. The TLSs generate random telegraph noise 
(RTN), which with a wide distribution of switching rates, 
gives 1/f nois d 22 * 2 ^ and behave very much as classi- 
cal stochastic variables such as fluctuating Ising spins. 
A quasi-Hamiltonian formalism was introduced recently 
which was found to be extremely suitable and versatile 
for describing the non-unitary temporal evolution of a 
quantum system acted on by such a classical stochastic 
process ^"^H 

However, in a circuit QED type of a setup, the qubit 
is not only coupled to the TLSs, but also to a cavity 
mode. This further complicates the picture. In the 
Jaynes cummings model, the offdiagonal dipole coupling 
term between the superconducting qubit and cavity can 
be roughly expressed as A oc \i • B rra \ — \\i being the mag- 
netic moment of the qubit and B rms is the rms value of 
photon's magnetic field). Now if a single TLS has a ran- 
dom time dependent fluctuating magnetic moment mj (t) , 
then it can be argued that the dipole coupling term will 
now look like A ex [fi + Trij(i)] • B rms . Thus fluctuations 
due to pure dephasing noise will induce fluctuations in 
the dipole coupling term which then introduces an en- 
ergy relaxation mechanism. 

The experimentally observed relevant energy scales are 
close enough for these effects to be noticeable. In case of 
a phase qubit, the vacuum Rabi splitting energy is about 
A = 100 MHz^(upto a factor of ft). While t he estimated 
TLS's spliting energy is about 10 MHz 11 26 and has been 
reported to be even as high as 45 MHz^. 

In this paper, we extend a recently developed quasi 
Hamiltonian formalisrrP^HII] to treat a qubit-resonator 
system that is affected by the presence of multiple TLSs. 
This system schematically shown in Fig{l] In the quasi- 
Hamiltonian method, the non-unitary temporal evolution 
of a quantum system acted on by a classical stochas- 
tic process in described in terms of the evolution of a 
non-Hermitian Hamiltonian. Decay processes in open 
quantum systems have been treated using using non- 
Hermitian Hamiltonians in the pastpSHMl 

Quite commonly, dissipative open quantum systems 
are treated using the Lindblad master equation 35 . The 
Lindblad formalism can be extended to quantum op- 
tics type systemP^ 3 -^ and can be solved by the use of 
methods such the stochastic wavefunction method^EU 
when many degrees of freedom are wanted. The stochas- 
tic wavefunction method has become a popular numer- 
ical technique since in comparison to the reduced den- 
sity matrix method as it provides additional information 
about the state of the systeirM H owever, t his and other 
Monte-Carlo wave function method j 36 l 38 l 42 l are somewhat 
numerically intensive. 

In our formalism, the dissipative dynamics of the en- 
tire qubit-resonator-fluctuator system is obtained by a 



single shot calculation. At the same time, we are able 
to individually set noise parameters for each fluctua- 
te which enables us to introduce a broad distribution 
of relaxation rates (which lead to 1/f noise) and TLS 
splitting energies. Another advantage of our method is 
that temperature dependence and spin-spin interactions 
(such as RKKY for the fluctuators), can be easily intro- 
duced - both of which are known to affect 1/f noise in 

squidPeh. 

Or calculations shows that the system displays some 
unusual characteristics due to the stochastic fluctuations. 
For example, behavior that closely resembles collapse and 
revival of Rabi oscillations occur at times that depend 
on the noise coupling strength. A broad distribution of 
noise coupling strengths leads to a broad distribution of 
this type of behavior over time. It is also seen that, par- 
ticularly for non-Markovian noise, the quantum photon 
bath can either speed up or slow down the relaxation and 
dephasing processes depending on the initial conditions. 

This paper is organized as follows. In sec (H] we discuss 
the model and our calculation method using the quasi 
Hamiltonian formalism within the Jaynes Cumming s 



model. In sec |III| we discuss we present approximate an- 
alytic solutions for the case of a single TLS. The effect 
of mutiple TLS with fixed or random noise strength dis- 
tribution and with a 1/f noise power spectrum is then 
discussed. This is followed by a particular example of a 
phase qubit in sec jlll C[ Lastly we present our summary. 



II. MODEL AND METHODOLOGY 

Consider the schematic of a superconducting(SC) 
qubit coupled to a resonator , as shown in figJTJ At the 
metal insulator interface there exist a number of spins 
which randomly flip at different instances of time. The 
Hamiltonian describing the entire system is as follows 



H = H s + H h + Hi 



sb 



(1) 



The system Hamiltonian, H s , for a single qubit (or a 
two level atom) in the presence of a single mode quantized 
electromagnetic field is given by the Rabi model. In the 
rotating wave approximation, or the Jaynes- Cummings 
model, only the energy conserving interaction term is re- 
tained 



H s = ^ycr 2 +uj'<Ja+ ^(aa + +a t cr_) 



(2) 



where, uj /2 is the energy separation between the excited 
state |1), and the ground state, |0), of the qubit and is 
proportional to the applied magnetic field. Note that 
we have set ft = 1. Here uj' is the photon frequency, a) 
and a are the photon energy creation and annihilation 
operators for a single mode, a± = a x ± ia y . In its most 
general form, A ex (0|d|l) is the dipole coupling term 
between the qubit states and d is the dipole operator. 
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The exact form of A depends of the SC qubit type and 
resonator type. 

The fluctuating two level systems at the metal insula- 
tor interface are modeled as flip- flopping Ising spins (fluc- 
tuates). If the fuctuators are statistically independent, 
then the noise Hamiltonian can be expressed in terms 
of the sum of contributions from individual fluctuations 
acting on the qubit 



U) 



(3) 



where, Sj(t) is a two level sequence that switches between 
±1 at random intervals of time, g^') is the noise vector for 
the j th fluctuate and |g| = g is the noise strength. Their 
autocorrelation function is (si(ti)sj(t2)} oc exp(— 2jj \t\ — 
^D^ijj 7j is the switching rate of the j th fluctuating Ising 
spin and each fluctuate has a Lorentzian power spec- 
trum. 

For our calculations we assume that the SC qubit's 
working point is such that, by it self, the SC qubit only 
has pure dephasing noise, i.e. gi = [0,0,gf]. The noise 
Hamiltonian then reduced to Ht> = ^2si(t)gfa z . 

This noise Hamiltonian causes fluctuations in the en- 
ergy levels of the qubit, which perturbs the wave func- 
tion |1) — >• . This in turn will cause fluctuations in 
the dipole coupling term A ex (0|d|l)'. This fluctuating 
dipole term will result in the presence of an off-diagonal 
energy relaxation term in the 2x2 subspace of the Jaynes 
Cummings Hamiltonian. 



H *b = 7,^2 9x ] ) s j (t)(aa+ + a 1 'cr_) 



(4) 



where gi^ oc (0 1 1) r — A and has exactly the same time 



dependence as g^J\ An alternate form of g^i\ more suit- 



able for a phase/flux qubit, is discussed in sec. IIIC 



More explicitly the full Hamiltonian in the {|0;n 
1), |1; n)} subspace is 



E 



raj + 



n+1 



n + 1 noj + + A - 



(5) 



where A 



uj is the field detuning. 



A. Transfer Matrix for the Qubit and the Photons 



the initial density operators for the qubit and the photon 
field, respectively, are 

MO) = |0> <0| (6) 
MO) = |*F(0)><* F (0)| = £|C„| 2 |n>H (7) 

n 

The time dependent density matrix is 

p(dt) = U[p Q (0)®p F (0)]U^ (8) 

where U = ex.p(—iHdt). We next take a partial trace 
over the field to obtain the reduced density operator of 
the qubit 

p Q (dt) = Tr F [p(dt)} (9) 
= ^(n\U\p Q (0)®p F (0)]Ul\n). (10) 



In general, the density operator for a single qubit can 
be written as 



p(t) 



1+ ^2 Vk(t)er k 

k=x,y,z 



(ii) 



where, rjk(t) are components of the Bloch vector, rj(t), 
and 1 77 1 is a measure of purity. 

The temporal dynamics of the quantum system can be 
reformulated as a transfer matrix equation by using the 
identity, given by EqJTTJ on both sides of Eq. [l0| 

I + 5> fc (dt)<7 fc = ^(n|£/(/(g)p F (0)£/t|n) (12) 

k n 

+ J2( n \ U lVk(0)(T k ® p F (0)]tf\n). 

n,k 

Noting that Tr[p F (0)} = J2 \ C n\ 2 = 1, we obtain 

J2vk(dt)a k = ^(nplrikiOfck® PF(0)]U*\n) (13) 

k n,k 

= Y,\C n \ 2 rik(P)U n (T k Ui (14) 



n.k 



Also note that the unitary matrix, J7, is block diagonal, 



where each of its 2 x 2 blocks are U n = e 



. Multi- 



plying both sides of Eq(l4] by aj and using the identity 
Tr[ak(Jj] — 25kj, the following transfer matrix equation 
can be obtained 



As we have mixed states, we seek solutions in terms 
of the density operator. We assume that at t = 0, the 
qubit and the single cavity mode is not entangled and 
are pure states, i.e. the system's initial density matrix 
is p(0) = Pq(0) (8) Pf(0). We can also further assume 
that the qubits are initially in the |0) state. Therefore, 



rj(dt) = T • 77(0) (15) 
where the elements of the 3x3 transfer matrix T are 

T kj = \ C n?Tr [UnCTkUlaj] (16) 

n 



4 



As an example, in the presence of only a single TLS(with 
s(t) = ±1), the two possible transfer matrices in the 
small time approximation are 



T±(dt) = ^\C n \ 2 {I + L z (ATg z )dt 



(17) 



L X (X ± g x )Vn + ldt} 



Similarly, for N TLS, there are 2^ possible transfer ma- 
trices T[ ±5±v>>± ]. 



Quasi Hamiltonian in the Presence of Many 
TLS 



The overall temporal evolution of the quantum system 
(single cavity mode plus qubits) is governed by the time 
ordered product of the transfer matrices obtained in the 
previous section 



r)(t) 



t/dt 



J7(0). 



(18) 



Whereas the classical stochastic noise process, for the 
two level fluctuators, is governed by the master equa- 
tion, W(t) = VW(t) 43 , where V is a matrix of transi- 
tion rates (such that the sum of each of its columns is 
zero) and W is the flipping probability matrix for the 
two level systems. For N uncorrelated fluctuators, the 
flipping probability matrix is: 



W = Wi<g>W 2 <8>...W 



N 



(19) 



where the flipping probability matrix for a single fluctu- 
ate, assuming equal occupation probability of the two 

states, is 



Wj(t) = 



-2 1} t 

1 - e" 2 ^'* 1 + e" 2 ^* 



(20) 



Here jj is the rate at which the fluctuator sj(t) switches 
between ±1. 

In a manner similar to path integrals, one can therefore 
construct a combined transfer matrix that describes the 
small time(dt) evolution of the quantum system coupled 
to the stochastic system (from the two types of transfer 
matrices discussed here) as follows 



W T. 



(21) 



Where denotes a Hadamard product and T is a square 
matrix (n being the number of qubits), each of whose 
columns consists of the lexicographically ordered trans- 
fer matrices [T[ +j+} . i+ ],T[ +}+} „_]...T[_ j _ } .._]] (the sub- 
scripts denote the TLS configuration). 

The state of the quantum system along with N clas- 
sical TLSs ,at time t, is given by T m = r m r m _i...ri, 
where m = t/dt. By equating T m = exp(—iH q t), one 



can obtain H n 



lim i(T 

dt^O 



I)/dt in the small time limit. 



Where H q is a time-independent non-Hermitian quasi- 
H audit oniarPSl. Its non Hermiticity implies that the time 
evolution operator ex.-p(—iH q t) need not be unitary, and 
hence can be used to treat dissipation in open quantum 
systems. 

Finally we arrive at our main set of tools for calculating 
dissipation in this system. In the presence of TV fluctuat- 
ing TLSs and a single cavity mode, the quasi Hamiltonian 
is 



±± q 



where 



Hq"f Hqg ~\~ Hq^j 



N 



3 = 1 
N 



(22) 



(23) 



H qg = i J2 T * J) ® (di j) L z + g^V^+lL x ) (24) 
i=i 

(25) 



H, 



quo 



il (AL Z - Wn + IL X ). 



Here, Li =XiViZ are the SO (3) -generators, Lq is a 3 x 3 
identity matrix, / is a 2^ x 2^ identity matrix and it is 
implied that 



Si) 



(i+i) 



» 



(26) 



Finally, the time dependent Bloch vector for the qubit 
is obtained using the following projection 

V(t) = (/iv|...(/2|®(/i|^|C' n | 2 e-^ n) *|H)®K 2 )...|^)7 7 (0) 

n 

(27) 

where \i) and |/) are the initial and final state vectors 
for a single TLS, that satisfy W\i(f)) = \i(f)) (these 
correspond to the zero eigenvalue solution of V). For an 
unbiased TLS (i.e. with equal occupation probabilities), 
|z) = |/) = [l,l]/V / 2. 



III. RESULTS AND DISCUSSION 

A. Single Fluctuator 

We begin our discussion by considering a qubit coupled 
to a resonator under the influence of a single fluctuator. 
For a single cavity mode initially in a coherent state, 



\C n \ 2 =exp(-(n»- 



(28) 



where (n) is the average photon number and \C n \ 2 is the 
probability that there are n photons present at t = 0. 

In order to obtain a clearer understanding of the Bloch 
vector's dissipative dynamics in such an environment, 
we solve Eq(27| numerically and obtain analytical expres- 
sions as well. The matrix exponential in the time evolu- 
tion operator can be approximated using the Zassenhaus 
expansion^ as follows 



e -iH q t _ e -iH' q t e -iH qu t e [H' q ,H qu ]t 2 /2 



(29) 
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where H' = H qi + H qg and the commutator [H' , H qcj ] = 
— (g z X + g x A)a z 0L r If we further assume that the 
resonator and the qubit are in resonance (A = 0), and 
noting that Eq 29 is only valid at short times, t, or when 
g z X is small. 

We obtain approximate analytic expressions for the 
Bloch vector components for two different sets of initial 
conditions. If the Bloch vector initially pointing in the 
x-direction, (i.e. rj Q = [1,0,0]), then the following com- 
ponents are obtained 



oo 2 

J2\Cn\ 2 [cos(\ n g z t 2 ) 9 "' 

n=0 



9xn + 9z 



(30) 



+g z sm(X n g z t 2 ) sin(2A n t) 



sin(Ot) 



-it 



n 



= ^2 \C n \ 2 g xn s'm(X n g z t 2 )cos(2X n t) 



sin(fit) 



-7* 



(31) 



n=0 



^2\C n \ 2 gxng z cos(X n g z t 2 ) 



2^ i - m 



9xn + 9 2 z 



(32) 



yxn 

sm(X n g z t 2 ) sm(2X n t) 



sin(Qt) 



- 7 t 



where 



C(t) = cos(fit) + ^ sin(fit) 



-7* 



(33) 



and n Vdxn+9 2 ~ 7 2 , flkn = ^V^Tn and A n 

Whereas if we start with r] = [0, 0, 1], then the follow- 
ing components are obtained 



Vx = X/ \ c n\ 2 9xn9z cos(2A n t) cos(X n g z t 2 ) 



n=0 



(34) 



In contrast if the Bloch vector is initialized along the x- 
axis (f] = [1,0,0]) then it will undergo either an oscilla- 
tory type decay or a monotonic decay to the center to the 
Bloch sphere as shown Fig(2]-a. From Eqs 3T|f33| it appar- 
ent that rj x (t) i s the only non-zero term if g xn = X n = 0. 

In Fig(2^i, the decoherence of the Bloch vector is shown 
in three different noise coupling regimes - in the strong 
coupling limit (non-Markovian noise, g z > 7), in the in- 
termediate regime (g z ~ 7) and in the weak coupling 
limit (Markovian noise g z < 7). \rj\ is a measure of the 
purity of the state and is zero only at the center of the 
Bloch sphere. 

As one crosses over from the non-Markovian to the 
Markovian noise regime, the Bloch vector's decay goes 
from being oscillatory to be monotonic because the 



trigonometric functions in ((t) (see Eq 33) become hy- 
perbolic functions for Markovian noise (7 > g z + g X n)- 
For non-Markovian noise, random Bloch vectors oscil- 
lations are caused by the g z s(t)a z term. For a given 
noise realization, the randomly switching s(t) causes ro- 
tations about the z-axis, however, as ensemble averaging 
restores chiral symmetry, the averaged Bloch vector al- 
ways travels in a straight line from pole to pole on the 
Bloch sphere^ (where the poles are |+) = (|1) + \0))/y/2 
and |— ) = (|1) — \0))/y/2) before diminishing to the cen- 
ter. 

Next, consider the case where the qubit is strongly cou- 
pled to the resonator and is completely decoupled from 
the fluctuator. The single cavity mode is initially pre- 
pared in a coherent state. In this case if r] = [1,0,0], 
then the Bloch vector will remain oriented along the 
x-direction indefinitely (which is again verifiable from 
Eqs 31 [33). Whereas if the Bloch vector is initialized 
along the z-direction then it will undergo the well known 
phenomenon of collapse and revival (CR) of Rabi oscil- 
lations as shown in Figj2}}. This CR process continues 
indefinitely with each revival being smaller in amplitude 



< = ^|C„| 2 [sin(2A n t)cos(A n ^ 2 )C(*) ( 35 ) 



+g z on^gj^—^-Lr-V 



ft 



Vz 



OO 

£ 

n=0 



C n \ cos(2X n t)cos(X n g z t ) 



2\9z 



9 2 xnC(t) 



9x 



■9i 



(36) 



Next, consider the Bloch vector's dissipative dynam- 
ics for two extreme cases. Note, that the Bloch vector 
dynamics shown in all the figures are obtained by nu- 
merically solving Eq|27[ 

In the first case let the qubit be subjected only to pure 
dephasing random telegraphic noise from a single fluc- 
tuator with negligible coupling to the resonator. Under 
these conditions, if the Bloch vector is initialized along 
the z-direction, then it will remain there indefinitely (as 
, where rj f z (t) 



seen in Eqs. 

if gxn ~ 



34 



36 

or 



land = 



77 =[1,0,0], y = 0.1, A, = 10"' 



77 =[0,0,1], g = 0,^ = 1 




x 10 



FIG. 2. Two extreme cases of the Bloch vector's dissipative 
dynamics, (a) TLS induced decoherence in the Markovian 
(7 > 9z\ intermediate and non-Markovian noise regimes with 
negligible coupling to the resonator, (b) Collapse and revival 
type behavior of Rabi oscillations induced by (n) = 10 pho- 
tons in the coherent state and in the absence of coupling to 
the TLS. Note the different initial conditions (r] ) in each case. 
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and less distinct from the preceding collapse. Increas- 
ing (n) results in more rapid collapses and more distinct 
revivals that are spaced further apart in time. 

These behavior can be explained as follows. In the 
limit of vanishing g z and g xn , Eqj36] reduces to the well 
known form: rj' z (t) = ^2 \C n \ 2 cos(2A n t). Each term 
in the summation (over n) represents Rabi oscillations 
(weighted by |C n | 2 ) that are associated with a definite 
value of n. At t = these different Rabi oscillation terms 
are all correlated. As time increases the destructive in- 
terference between these weighted oscillatory terms leads 
to the collapse of Rabi oscillations. However, at longer 
times the \C n \ 2 cos(2A n t) type terms will be in phase 
again and their constructive interference leads the the 
revival of Rabi oscillations. The revival phenomenon is 
of course is a purely quantum mechanical feature and is 
dependent on the photon distribution C n . If one were to 
take a continuous distribution of photons, then no revival 
occurs after the collapse and the quantum dissipative pro- 
cess then becomes indistinguishable from that due to a 
classical stochastic process. 

Next, we consider the case where the qubit is subjected 
to a single fluctuator and is also coupled to the resonator 
with (n) = 10. As discussed in secjnj in the presence of 
a pure dephasing noise, the fluctuating dipole coupling 
term(A) induces energy relaxation with noise strength 
g xn . For our calculations we take g xn = g z ^n\jX n + -0- 

The presence of g xn makes the picture more compli- 
cated as seen in Figj3j-a, where the decoherence of the 
Bloch vector is shown in non-Mar kovian, intermediate 
and Markovian noise regimes for r] = [1,0,0]. It is 
seen that the decoherence process is significantly slowed 
down for Markovian noise and noise in the intermediate 
coupling regime. Whereas for non-Mar kovian noise, this 
slowing down of the decoherence is less significant. When 
the initial conditions are changed to r] = [1,0,0], then 
there is no slowing down of decoherence as seen in Fig(3[ 
b. On the contrary the coherence times terms deteriorate 
with increasing noise coupling strength, g z . However, 



certain unusual collapse and revival type phenomenon oc- 
curs in the non-Mar kovian noise regime. While typically 
the occurrence of revivals (as a function of At) depends 
only on the photon distribution |C n | 2 , in the presence of 
a fluctuator, smaller revivals can occur at earlier times 
depending on g z \ As explained earlier, typically con- 
structive interference between cos(2A n t) type terms (as 
permitted by C n ) leads to revivals for a pure resonator- 
qubit system. However the presence of the cos (g z At 2 ) 
type terms in the expansion of the Bloch vector compo- 
nents (Eqs. (34j36|, will cause the constructive interfer- 
ence and hence the revival type phenomenon to occur 
earlier and with greater prominence, if g z is large (or for 
non-Mar kovian noise). 

An experimental observation of this phenomenon 
would be a good indication that TLS are present. If 
there are only a small number of TLS, then then a fit 
of the theory of this section to the results can reveal the 
coupling strengths and switching rates of the TLS. 



B. Multiple Fluctuators 

Typically, most experimental systems have many TLS 
with a broad distribution of switching rates. In this sec- 
tion we extend our calculations to this case. 

For the following set of calculations we use a set of 
eight fluctuators (as Hq 1 ^ scales as 2 n , calculations with a 
large number of fluctuators quickly becomes prohibitive) . 
Each individual fluctuator has a Lorentzian power spec- 
trum and a broad distribution of their relaxation rates, 
7 (ranging from 0.14 to 2 x 10 -5 ), results in a 1// noise 
power spectrum as shown in Fig|4}a. For simplicity and 
to limit the number of parameters used (in this section) , 
we have held gjS the same for all the fluctuators. 

The decoherence of the Bloch vector for different noise 
strengths is shown in Fig. [4}b and c for two different sets 
of initial conditions, r} . As expected the overall rate of 
decoherence increases with increasing noise strength, gj. 



77 o =[1 ,0,0], y = 0.1 , X = 0.1 77 o =[0,0,1], y = 0.1 ,X = 1 




1 2 10 20 



FIG. 3. The Bloch vector's decoherence in the Markovian 
(7 > g z ) , intermediate and non-Markovian noise regimes with 
(n) = 10 photons in the coherent state for the Bloch vector 
initially in (a) j] Q = [1, 0, 0] with A = 0.1 and (b) r\ = [1, 0, 0], 
A = 1. These calculations are carried out numerically. 
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FIG. 4. (a) Power spectrum showing 1/f noise for a set of 
eight TLSs. (b) Decoherence of the Bloch vector for different 
noise strengths with A = 0.1, (n) = 10 coherent state photons 
for an initial state r] = [1, 0, 0] and (c) r\ — [0, 0, 1]. 
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FIG. 5. (a) Decoherence of the Bloch vector for different noise 
strengths with A = 1, (n) = 10 coherent state photons for an 
initial state r\ Q — [1,0,0] and (b) r\ — [0,0,1]. Note that 
the power spectrum for these fluctuators is the same as that 
shown in Fig|4|a. 



that on increasing A, strong suppression of decoherence 
is seen in the Markovian and intermediate noise regimes 
for T] = r\ x (see Figj5]-a). For non- Markovian noise the 
decoherence is more apparent, but still slower than the 
previous case of A = 0.1. For r] = rj Z: however, the dis- 
sipative behavior of the Bloch vector shows much richer 
behavior (see Figj5]-b). The unusual collapse and revival 
type phenomenon seen earlier for the single fluctuate 
case (see Figj3]-b) occurs here as well for gj = 1. How- 
ever for this to be more apparent in the case of multiple 
fluctuators, the noise must be strongly non-Markovian 
and the qubit should couple to the resonator strongly. 
Overall, A should be large enough to overcome the sup- 
pression of this collapse and revival type phenomenon by 



C. A Phase Qubit Example 



The variations in the dissipation rate for different gj is 
most apparent if one starts with an initial state of rj x as 
shown in Fig. [4[b. For a full spectrum of fluctuators, if 
gj is smaller than the smallest 7 then the noise due to the 
TLSs falls purely in the Markovian noise regime, where 
the rate of dissipation will be the least. If gj falls some- 
where in between the selected range of 7s, then one has a 
mixture of Markovian and non-Markovian noise sources, 
which is the case for gj = 0.01. However in this case, 
the Markovian noise sources tend to dominate and the 
oscillatory behavior (as apparent from Eq|33|), typically 
seen for non-Markovian noise, tends to get washed out, 
as shown for an initial state of r] z (Fig. |4j-c). However 
now, collapse and revival of Rabi oscillations are seen at 
longer times due to the presence of a single photon mode, 
similar to that of Figj3|b (this is not shown in Fig. [4j as 
we wish to focus on the initial collapse). 

In the case of gj =0.1 one has a mixture of Markovian 
and mostly intermediate noise sources. This leads to an 
oscillatory type of behavior where small oscillations are 
superposed on top of a smoothly decaying function. Now, 
if g 3 ; = 1, then one is entirely in the non-Markovian noise 
regime, and the decay of the Bloch vector is strongly 
oscillatory. However, now the decay of the Bloch vec- 
tor is far more strongly affected by the resonator due to 
cos(g z \t 2 ) type terms in Eqs. 31p6 The initial collapse 
of Rabi oscillations now mixes with TLSs induced oscil- 
lations with characteristic frequencies of flj. A broad ^ 
distribution of Qj (due to 7^) results in the complicated 
beating of the Bloch vector seen in FigE] (for gj = 1). 
However now, the subsequent revival of Rabi oscillations 
(similar to that of Figj3 b for gj = 1) will not be seen 
at longer times as this behavior will be suppressed by 
e -J2ijt due to the presence of multiple fluctuators. For 
this to be visible we have to increase the resonator cou- 
pling strength, which is what is done for the next set of 
calculations. 

FigjF] shows the dissipation of the Bloch vector for 
A = 1 for two sets of initial conditions. It is apparent 



Our discussion of the TLS induced Bloch vector dy- 
namics, in the previous section, was general and did not 
pertain to any particular type of qubit. In this section 
we consider the example of a phase qubit. In general, all 
superconducting qubits are characterized by the nonlin- 
earity of its Josephson junctions, which can be thought 
of as a nonlinear inductor. When coupled to a capacitor 
(either external or arising from the junction), it forms a 
nonlinear LC oscillator. In the case of phase qubits, it so 
happens that the characteristic impedance of this non- 
linear LC oscillator matches that of microwave transmis- 
sion lines/resonators, which then allows for strong qubit- 
photon coupling. 

The general Hamiltonian of a currant biased phase 
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FIG. 6. Decoherence of the Bloch vector for A = 100 
MHz, g xo = 10 -2 , (n) = 10, r\ — [0, 0, 1], and a normal distri- 
bution of gj with mean (g) = 10 MHz and standard deviation 
of (a) s = 0.001 MHz, (b) q = 0.01 MHz and (c) q = 0.1 MHz. 
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FIG. 7. Mean value, (g z ), as a function of the location of the 
peak values in the fourier spectrum (see Fig[6]d-f) frequency 
for various (n). These calculations were carried out for seven 
TLS with <^ = 0.1 MHz and k = 0.01. 



qubit is H<p = q 2 /2C - I $ (cos(p - I b p/I )/2ir. Where 
q is the charge operator, C is the capacitance, <£ is 
the magnetic flux quantum, I b is the bias currant of the 
Josephson junction(JJ), I Q is its critical value and cp is the 
phase across the J J. By suitably biasing the J J one ob- 
tains a washboard type potential energy landscape where 
the qubit states correspond to the two lowest quantized 
energy levels. This energy difference corresponds to uj 
m Eq|l] In our calculations, however, this only appears 
in the form of the level detuning A, when calculating the 
dissipative dynamics (see Eq .25). For all of our calcula- 
tions we take A = 0. And for the vacuum Rabi split- 
ting energy, we use the experimental value of A = 100 
MHz 28 (note that H= 1). 

Experimentally it is seen that the TLS have a gaus- 
sian distribution of splitting energies with a relatively 
small variance 29 for most fluctuates. Their estimated 
mean fluctuator splitting energies were about 10 MHz, 
however a few of the gjs were as high as about 40 MHz. 
Other estimates of gj obtained by fit ting to spin echo 
data vary between about 9.6 to values as high 

as 135 MH d 26 * 4 ^ . Though the various experimental esti- 
mates vary, it is quite clear that the TLS splitting energy 
is much higher than a factor of gL^B (where gL is the 
Lande ^-factor and he is the Bhor magneton) for an elec- 
tron spin. 

In the previous section gj was assumed to be the same 
for all the fluctuates. Here we examine how the collapse 
and revival type behavior would be affected if one were 
to consider a number of fluctuates with a random distri- 
bution of gj . For the next set of calculations we consider 
a set of six Gaussian distributed random TLSs with a 
mean value of (g z ) = 10 MHz and a standard deviation 
of = 0.01,0.1 and 1 MHz. 

As mentioned earlier, the vacuum Rabi frequency A oc 
H • B rm p^{ where fi is the magnetic moment of the qubit 
and B rms is the root-mean-square value of photon's mag- 
netic field). In the presence of multiple TLS with mag- 
netic moment m^, we argue that the dipole coupling 



term will look like A ex [/i + J2j m j^{^)\ ' Brms- From 
this, the fluctuating coupling term can be written as 
gin = Kgi^y/i + n, where k is phenomenological scal- 
ing constant and we set n = 0.01. Now, note that all of 
our calculations, the collapse and revival type behavior 
is invariant with respect to the Kt time scale. Hence the 
noise induced Bloch vector dynamics will be observable 
even for the smallest k. 

The resultant temporal dynamics and the respective 
Fourier transforms are shown in FigjHJ The first peak 
in the Fourier spectrum corresponds to the initial set of 
oscillations while the second less prominent peak is due 
to the revived secondary set of oscillations. The width 
of these peaks corresponds to the various frequency com- 
ponents of these oscillations. As seen in the figure, the 
collapse and revival type phenomenon is more apparent 
if ^ is small. For a larger ^, the visibility of this effect 
diminishes due to the superposition of the collapse and 
revival type oscillations with more widely varying fre- 
quencies. This washing out effect is also apparent in the 
Fourier spectra, particularly in the significant broaden- 
ing and lowering of the secondary peak. However, we also 
see that this washing out effect ,with increasing can be 
countered to some extent by increasing (n). 

The dependency of the frequency of these CR oscilla- 
tions on (g z ) and on (n) is examined more closely nest. 
In FigjTJ the frequency at which the peak in the Fourier 
spectra occurs, uo pe ak ? (corresponding to the first set 
of oscillations) is shown as a function of (g z ) and (n). 
As expected uj pea k varies linearly with (gj) and as (n) 
is increased uj pea kS dependence on (g z ) becomes more 
discernable. 

In general the dependence of uj pea k on (g z ) can be 
somewhat complicated. However, based on our analytic 
solutions(for e.g. see Eq{36|), we found the following re- 
lation to be approximately true. 



^peak 



(37) 



Even more encouragingly, this relation is independent 
of the number of TLSs (which was explicitly verified 
for upto seven particles). Furthermore the accuracy of 
Eq|37| improves with increasing average number of co- 
herent photons (n). 

This has important consequences from an experimen- 
tal point of view. This implies that the mean value of 
the TLS splitting energy could be determined by load- 
ing different number of photons into the resonator and 
observing the frequency shift in 0J pea k- The only detri- 
mental factor to this approach would be a vary large 
standard deviation, in the distribution of #jS, which 
would make the Fourier spectrum less discernable. This 
however could still lead to further physical insight into 
the system. 
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IV. SUMMARY 

In summary we have calculated the dissipative dynam- 
ics for a superconducting qubit coupled to a resonator, 
where the qubit is affected by a number of fluctuating 
two level systems with a 1/f noise power spectrum. The 
decoherence of the Bloch vector is calculated using a 
quasi Hamiltonian formalism whithin the Jaynes Cum- 
mings model. 

We argue that fluctuations that cause pure dephasing 
noise will also induce an energy relaxation mechanism 
once the qubit is coupled to the resonator. The simul- 
taneous coupling of a qubit to a coherent photon bath 
and classical TLS is shown to result unusual dissipative 
behavior. It is seen that collapse and revival of Rabi 
oscillations occur at times that depend on the noise cou- 
pling strength. This collapse and revival type behavior 
is particularly strong in the non-Markovian noise regime. 



Depending on the initial orientation of the Bloch vector, 
the coherent quantum photon bath either speeds up or 
slows down decoherence, which is also particularly strong 
for non-Markovian noise. 

A broad distribution of noise coupling strengths leads 
to the broadening of the collapse and revival type behav- 
ior, which makes it less apparent. The frequency of these 
Rabi oscillations depends on the mean noise strength. 
We find an approximate analytic relation for this find it 
to be independent of the number of fluctuators. From 
an experimental point of view, it is perceivable that with 
controllable A and (n) one can estimate the average noise 
strength and various characteristics of its frequency spec- 
trum. 

The authors wish to thank Robert McDermott for a 
number of useful discussions. This work is supported by 
DARPA-QuEst Grant No. MSN1 18850. 
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